Extremely Sub-wavelength Planar Magnetic Metamaterials 
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We present highly sub- wavelength magnetic metamaterials designed for operation at radio frequen- 
cies (RFs) . A dual layer design consisting of independent planar spiral elements enables experimental 
demonstration of a unit cell size (a) that is ~700 times smaller than the resonant wavelength (Ao). 
Simulations indicate that utilization of a conductive via to connect spiral layers permits further op- 
timization and we achieve a unit cell that is Xo/a ~ 2000. Magnetic metamaterials are characterized 
by a novel time domain method which permits determination of the complex magnetic response. 
Numerical simulations are performed to support experimental data and we find excellent agreement. 
These new designs make metamaterial low frequency experimental investigations practical and sug- 
gest their use for study of magneto-inductive waves, levitation, and further enable potential RF 
applications. 



Experimental verification of a negative index of 
refraction [ll Q confirmed a four decade old theoretical 
prediction [3] and initiated a new area of research in arti- 
ficial electromagnetic (EM) materials known as meta- 
materials (MMs). In the last decade metamaterials re- 
search has yielded many exotic effects including: in- 
visibility cloaking^, perfect lensing^, 0], and perfect 
absorbers 0. Interest in metamaterials stems from their 
available use in nearly any band of the electromagnetic 
spectrum and their ability to achieve almost any desired 
response [§-11 1 . 

While there are relatively few technical restrictions 
limiting operational frequencies of metamaterials there 
are, however, some practical hurdles. To date there has 
been a noticeable lack of metamaterials demonstrated at 
low frequencies - on the order of 100 MHz or lower. Al- 
though metamaterials may still be of use in this frequency 
range, the main drawback is that the dimensions can be- 
come impractically large. Wavelengths for radio frequen- 
cies range from thousands of meters (for 100 's kHz) to 
lengths of meters (100's MHz). Since it is not uncom- 
mon for individual metamaterial elements to have sizes 
on the order of ~ Ao/10 — Ao/100, (where Ao is free space 
wavelength), the side length of a single RF metamate- 
rial element could be expected be around 100 m to 10 
cm. In many situations this size may be unacceptable 
for real-world applications, particularly so for 'portable' 
devices. The development of a new class of metamateri- 
als with a significant reduction in the ratio of its phys- 
ical size 'a' to its resonant wavelength Ao may create 
new opportunities for low frequency metamaterial-based 
applications 12 - 15| . 

Magnetic properties exhibited by natural materials are 
typically weaker and less common than their electric 
counterpart. This is due to the fundamental difference in 
how magnetic and electric responses are generated - the 
former being produced either from orbital currents or in- 
trinsic unpaired spins. This fundamental lack of diversity 
in materials hinders full development of electromagnetic 
devices. Metamaterials with a //-negative response have 
been shown to improve the effectiveness of systems that 



make use of magnetic fields at radio frequencies such as 
increasing the resolution of MR jJla , Jl7| and enhancing 
wireless power transfer efficiency [18|, ll9(- 

Here we propose that dual layer metamaterial spi- 
rals may be used in order to achieve extremely sub- 
wavelength RF magnetic metamaterials. We utilize and 
detail a time domain method for characterization of the 
complex magnetic response (h(lj)) of the RF metamate- 
rials and support experimental results with simulation. 

As a starting point for exploration of RF mag- 
netic metamaterials we begin with a design based on 
a planar s pira l structure capable of achieving large 
inductances 2(3- 2^ - as it is trivial to wind a significant 
length of wire into a small area. The spiral response 
can be approximated based on a simple LC resonator 
model uiq = l/V LC, where ujq is the resonant angular 
frequency, L is the inductance, and C is the capacitance 
of the MM. Additional reduction of the unit cell size 
may thus be achieved by maximizing L and C simul- 
taneously. In single layer planar spirals the capacitance 
arises solely due to interactions between adjacent metallic 
windings 23]. However, since typical lithographic meth- 
ods produce extremely thin metallic layers the capaci- 
tive values obtained in this manner are small. We may 
drastically increase the capacitance by adding a second 
spiral layer, thus creating an interlayer interaction be- 
tween the two planar structures. In this way both the 
larger area of the broadside metallic surfaces and a sup- 
porting dielectric in-between will provide an increased 
capacitance. A similar idea had been demonstrated in 
the microwave range, "the broadside-coupled" split ring 
resonators (SRRs) were found to have a resonance oc- 
curring at much lower frequency than the original single 
layer designj^Q. The addition of the second layer will 
also provide an enhancement of the metamaterial induc- 
tance if added in series. 

We investigate two different dual layer planar spirals 
- circular and square windings, see Fig. [1] (b). Simula- 
tions were performed with a commercially available 3D 
full wave finite element frequency domain solver. Excita- 
tion from a waveguide port is used where the electric field 
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FIG. 1: (Color online) (a) Schematic (exploded view) of the 
dual-layer square spiral metamaterial (top) and a side view 
showing material layers (bottom), (b) Schematic drawings of 
circular (top) and square spirals (bottom), (c) Schematic of 
the experimental apparatus as described in the text. 



is along the x direction, and the magnetic field is along z 
- see Fig. 1(a). The simulated circular spiral design has: 
21 turns, a linewidth w = 170 /im, line spacing g = 320 
/im, inner radius r^„ = 1.2 mm, outer radius r ou t =11.7 
mm, and unit cell size (a) of 25.4 mm. The square spiral 
has: 13 turns, a linewidth w = 210 /im, line spacing g = 
320 /im, side length of the central capacitive pad sq = 1.9 
mm, and the side length of the square coil is 25 mm with 
a unit cell size (a) of 26 mm. The dielectric slab for both 
patterns is 203 /im thick and was modeled using a dielec- 
tric value of 3.6 + 0.13i. All metallic components were 
modeled as copper with a conductivity of a = 5.8 x 10 7 
S/m. For all cases, a vacuum box (a x a x 40(mm)) 
was used, and the geometry of the spiral, including the 
size of the centrally located capacitive pad, was modified 
in order to obtain the greatest oscillator strengths. The 
simulated scattering parameters, S%i and S21, allow one 
to extract effecti ve p arameters via the transfer matrix 
method approach [U Ea]. 

The extracted permeability for each magnetic meta- 
material is shown as red dash curves in Fig. [2] (c,e) for 
the circular spirals and (d,f) for the square spirals. A 
prominent magnetic Lorentzian type resonance can be 
observed for each configuration described by, 
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where fi^ is the relative permeability at frequencies 



higher than the resonance, 



^ nip 



is the magnetic plasma 



frequency, ujq is the center frequency of the resonance, 
and 7 is the loss. 



In order to characterize the complex permeability 
(p(to) = /ii + 1/12) response of our magnetic metama- 
terials, (/ii and \Xi are the real and imaginary portions 
of the magnetic response function, respectively), we uti- 
lized a time-domain spectroscopic method. A schematic 
illustrating our experimental apparatus is shown in Fig. 
[T](c). In our setup a coil is driven by a voltage pulse from 
a function generator creating a transient magnetic field. 
A second coil (pickup coil) is co-axially aligned with the 
drive coil, (spaced a distance z away) and measures the 
induced time dependent electromotive force (em fit)) due 
to the incident magnetic field created by the drive coil. 
Based on Faraday's law of induction, the induced emf of 
the reference and the samples can be expressed as, 



emf{t) r = -Afj, ^- 
emf(t) s = -A/irMo^ 



(2) 



where A is the area of the pickup coil, /io is the free space 
permeability, and H z is the z component of H field. 

The effective permeability (fi r ) of the samples can then 
be trivially derived from Eq. [5J when the areas A of the 
pickup coil and metamaterial samples are made equal. 



F{emf{t) s ) 
F(emf(t) r ) 



(3) 



In Eq. (3) F stands for the complex Fourier transform. 

The statement that the effective permeability can be 
described by Eq. 3 above carries with it the assumption 
that the time dependence of the H fields between the 
sample and reference cancel. Although this is a reason- 
able supposition we next show that this follows directly 
from Maxwell's equations and is verified by the experi- 
mental results. 

Ampere's law specifies that the bound current, which 
is related to the magnetization via Jb = V x M, and 
thus to the magnetic permeability, is what changes upon 
placing a magnetic material in the loop of the pickup coil. 
As H is only related to the true free current in Maxwell's 
Equations, this remains constant in our measurements 
and thus cancels for both the sample and reference mea- 
surements. Of course the above assumes that there is 
no change in the mutual coupling between the drive and 
pickup coils due to placement of the measured material, 
which would change the free current. We use the time 
domain method described here to characterize the com- 
plex magnetic response of metamaterials. In all cases we 
have made the areas of the pickup coil and metamaterial 
samples equal. 

Metamaterial samples were fabricated by using stan- 
dard printed circuit board (PCB) photolithography. All 
fabricated dimensions were identical to those simulated, 
and photographs are shown as the insets to Fig. [5J We 
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use the time domain method described above to charac- 
terize the complex magnetic response of metamaterials, 
see Fig. Q] (c). The time dependent magnetic field was 
generated using a coil antenna driven by a transient volt- 
age. The peak to peak voltage (V pp ) was set to 10 V for 
a period of 80ns with a duty cycle of 37.5%. A receiving 
coil spaced a distance of 8 cm away was connected to 
a digital sampling oscilloscope which captured the time 
dependent electromotive force induced by the incident 
magnetic flux. The shape and dimensions of the receiv- 
ing coil were chosen such that it matched that of a 3 x 
3 array of our magnetic metamaterials. For each sample 
we performed two measurements. The first was a refer- 
ence where nothing was placed inside the receive coil. For 
sample characterization we placed the metamaterial ar- 
ray into the center of the receiving coil. All measurement 
were collected with an average of 128 wave forms. 

Representative time pulses are shown in Fig. [2] (a) for 
both the sample (black curve) and reference measure- 
ment (red curve). A noticeable difference can be seen 
between the two measured time pulses as the reference 
signal has almost completely died out after 600ns, while 
the metamaterial signal 'rings' for a much longer time 
(over 1800ns). A plot of the frequency spectra calculated 
from the Fourier transform of the time data is shown in 
Fig. H] (b). Applying the extraction method of Eq. |3J 
outlined earlier, one can determine the effective perme- 
ability fi r (u>) of the magnetic metamaterials. 

The real and imaginary portions of the measured per- 
meability are plotted as grey curves in Fig. [5] for both 
the circular (top panels) and square (bottom panels) sam- 
ples. As can be observed both spiral metamaterial struc- 
tures yield magnetic resonances with Lorentz-like line- 
shapes. The circular spiral has its resonance centered 
at 17.3 MHz and the square spiral at 16.8 MHz. The 
real part of the measured relative permeability for both 
designs reaches a negative value of -0.25 for the circular 
design and -1.51 for the square design. The full width at 
half maximum (FWHM) is 0.554 MHz and 0.794 MHz 
for the circular and square designs, respectively. The 
calculated Q-factors are 30.3 and 21.8 and both designs 
achieve decent oscillator strength at resonance. The sim- 
ulated permeabilities (red dashed curves in Fig. [2][c)-(f)) 
are nearly identical to the measured ones. For example 
the simulated and experimental peak /i values differ by 
only 0.2% and 0.1% for the square and circular spirals, 
respectively. Some discrepancies in the amplitudes and 
center frequencies for both designs are observed. Never- 
theless, the experiments are in excellent agreement with 
the simulations. The errors can be attributed to imper- 
fections from the fabrication process. 

We compared experimental permeabilities to a 
frequency-dependent Drude-Lorentz oscillator, as shown 
in Eq. [TJ The black dotted curves of Fig. [2] are fits to 
the experimental data and we find good agreement with 
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FIG. 2: (Color online) (a) Time domain waveform recorded 
by the pickup coil for the metamaterial (black curve) and 
a reference signal (red curve), (b) The amplitude spectrum 
for the metamaterial (black) and reference (red) obtained by 
a Fourier transformation of (a). Experimental (grey curves), 
simulated (red dashed curves) and fits (black dotted curves) to 
the real (c,d) and imaginary (e,f) parts of [i r . Panels (c) and 
(e) show the real and imaginary permeability, respectively, 
for the circular spiral. Panels (d) and (f) show the real and 
imaginary permeability, respectively, for the square spiral. 



7 = 2tt x 0.745 (MHz), and w = 2tt x 17.3 (MHz) for 
the circular geometry, and /loo = 0.87, u mp = 2ir x 6.6 
(MHz), 7 = 2ttx 0.587 (MHz), and loq = 2ttx 16.8 (MHz) 
for the square geometry. 

Radio frequency magnetic metamaterials investigated 
here yield reasonably strong resonances although their 
physical size (a) is only a small fraction of their reso- 
nant wavelength. Circular spiral metamaterials achieved 
a = Ao/683 and square designs a = Ao/703. In com- 
parison to simulated results for a single layer spiral, (not 
shown), the resonant wavelengths of the dual layer de- 
signs are nearly six times smaller, although all have 
roughly the same unit cell size. It's worth noting that 
in-spite of the strong sub- wavelength sizes of metamateri- 
als presented here, there are simple improvements which 
will further increase Ao/a. For example we have chosen 
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FIG. 3: (Color online) (a) Simulated real (black) and imagi- 
nary (red) relative permeability for a dual layer circular spiral 
metamaterial with an interlayer via. (b) Experimental (solid 
curves) and simulated (dashed curves) relative permeability of 
a circular spiral metamaterial with conductive via - described 
in the text. 



FIG. 4: (Color online) Simulated results: (a) Measured per- 
meability for several values of coil separation z. The peak to 
valley amplitude (A pv ) of fi r (b) and the resonant frequency 
(ujo/2-k) (c) versus the coil separation distance (z). Solids red 
lines are linear fits to the data. 



to fabricate designs in which the two planar spirals are 
independent and only capacitively coupled through the 
substrate. We chose this geometry simply because these 
structures are easier to fabricate. However by connect- 
ing the two layers with a conductive via one may push 
sizes of these structures even more sub-wavelength, as 
the unwound length effectively doubles. 

Figure |31 presents simulations of the square spiral, with 
a connecting via, and a unit cell size a = 15.5 mm, 
linewidth w = 100 fj,m, line spacing g = 100 /im, board 
thickness of 203 /im, and consisting of 25 windings. Sim- 
ulations indicate this design should exhibit a strong mag- 
netic resonance at 9.34 MHz and the extracted complex 
fi is shown in Fig. [3] (a). For this geometry we find 
a = Ao/2072. Although the design discussed above is well 
within both laboratory and commercial fabricational ca- 
pabilities, here we chose to experimentally demonstrate a 
slightly less subwavelength structure, and simply modify 
the circular spiral shown in Fig. Q] (a) and (c). Holes 
were drilled through the center capacitive pad of each 
unit cell and a wire was soldered on both sides to create 
the connecting via. We then characterized the sample 
and extracted the complex permeability as shown in Fig. 
|3] (b). We find a resonant frequency of 8.94MHz and 
thus this modified spiral magnetic metamaterial yields 
a = Ao/1321. It should be noted that the experimen- 
tal structure presented here was by no means an opti- 
mized design, i.e. we did not fine tune capacitive and 
inductive reactances which could be responsible for the 
significantly lowered oscillator strength. Additionally a 
notable Ohmic loss resulted from the manner in which 
the conductive via was experimentally implemented. The 
metamaterial geometry was simulated where the via was 
modeled as a resistive element (R=13f2) and we find ex- 
cellent agreement with experiment as shown in Fig. [3] 
(b). 

As a verification to the validity of the experimental 
time-domain method presented here, we performed mea- 



surements of the response of magnetic metamaterials for 
several coil separation distances ranging from z = 4 cm 
to z = 24 cm. In all cases the sample under investigation 
was placed directly in the plane of the pick-up coil. As 
before all sample measurements were characterized and 
compared with a reference measurement. The induced 
maximum peak-to-peak value of the emf for the refer- 
ence measurements was observed to change from 2V for 
the closest spacing to 40 mV for the farthest. Calculated 
jj, r values are shown for three characteristic separations 
of z=8, 16, and 24 cm in Fig. H (a). In Fig. g] (b) and 
(c) we plot, respectively, the peak-to- valley amplitude of 
l^ r (u)) and the resonance frequency u>o, both as a func- 
tion of coil separation. The red horizontal lines of Fig. 
|4] (b,c) are fits to the data and have nearly zero slope 
to within 10 -4 . The measured center frequencies only 
deviate by a maximum of 1% of the average value for z 
distances studied here. Variation, however, in the peak- 
to- valley fi r values is as much as 5% though we attribute 
this discrepancy to noise in the system. 

We computationally and experimentally demonstrated 
extremely sub-wavelength planar magnetic metamateri- 
als. Experimental resonant wavelength to unit cell size 
values were shown to be as much as ^1300, and an ideal 
simulated design was shown to be Xo/a ~ 2000. We find 
excellent agreement between experimental and simulated 
magnetic permeabilities, for all structures investigated, 
and attain reasonable oscillator strengths. A novel char- 
acterization method has been presented capable of deter- 
mining the complex magnetic response of planar materi- 
als at RFs. Compact sub-wavelength designs presented 
here make low-frequency study of metamaterial practical 
which may lead to RF metamaterial applications. 
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